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Abstract 

We give a new proof of the well-known upper bound on the correlation function of a gas of 
non-overlapping dipoles of fixed length and discrete orientation working directly in the charge 
representation, instead of the more usual sine-Gordon representation. 
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1 Introduction 

+-> 

s 

Although dipole gases themselves do not exhibit the Kosterlitz-Thouiess transition, they have been 
used as a natural starting point for the study of this transition in Coulomb systems at low temper- 
^ ! atures and small activities Qj. 

Differently from the behavior of correlations in Coulomb systems, charge correlations in dipole 
gases have lower and upper bounds with power law decay. The lower bound is a consequence of 

_ _ i 

Jensen's inequality in the charge variables []J. The main ingredient to obtain an upper bound is 
the sine-Gordon representation of the gas. Using this representation Frohlich and Spencer have 

*Supported by FAPESP. Present Address: Department of Mathematics, University of California at Davis, Davis 
CA 95616-8633, USA 

■f Supported by CNPq. 



1 



proved a theorem: 

Theorem The two-point correlation function of external charges for dipoles with hard core is 
bounded above by 

G{x,y) <C p e- mloslx - yl (1.1) 

with C/3 < oo and 

m = -^-(l + ri + e)/^-^ 16 )- 1 
for any e > 0, where (3 is the inverse temperature and z is the activity. 

The correlation function above can be equivalently written as G(x,y) = e~@ Ve ff( x ' y \ where the 
effective potential V e ff(x,y) is the difference between the free energy of the system in the presence 
of the external charges and the free energy in the absence of the external charges. The physical 
content of the above theorem is that the long range dependence on the position of the external 
charges in fll.lj) implies that dipoles do not screen [Q. 

In this paper we dispense with the usual sine-Gordon representation and establish the upper bound 
fll.lj) for a gas of non-overlapping dipoles of fixed length and discrete orientation, working directly 
in the charge (or gas) representation. 

To obtain the above result we use an electrostatic inequality in the form of charge translations in 
order to drive down the activity of the dipoles and obtain estimates for small activities, followed by 
a direct cluster expansion of the partition function to extract the fall-off of the correlation. (The 
extension of the technique to Coulomb gases follows along the same lines). 

Our method provides an alternative procedure for dealing with systems for which no sine-Gordon 
representation is available. (See H for the use of a charge representation in connection with 
percolation models.) In particular, we hope that the charge representation will enable us to give a 
direct proof of the analyticity of the pressure both for dipole and Coulomb gases. 



2 Dipole Gases 

We consider a lattice dipole gas with discrete orientation in a finite region A. Let L = 4rZ 2 be a 
sublattice of A of mesh 4r, contained in a box A inside A and concentric to it. A configuration of a 
lattice dipole is described by a function JIl = {$(])}, where fl(j) = (/ii(j), /^(j)) is the total dipole 



moment at the site j G L. The dipoles are non-overlapping if their fixed length is r. For simplicity 

we shall take r — 1. 

The partition function is given by 

Z = £ e-^II^W)), (2.1) 

pL={P(j),jeL} j£L 

where the interaction energy is 

e (P-l) = Yl S ^{3)^{k)W ai {j, k) 

j,k a,y 

and the two-body potential is given by 

W ai (j, k) = d a d y C(j,k) 

where d a is the usual finite difference derivative in the a-direction, and C(j,k) is the Green's 
function of a lattice Coulomb potential in two dimensions. 
The "a priori" distribution d\ of is given by the measure 

d\(Jt) = 5(/I) + z5(\/I\ 2 - 1) 

where z is the activity, the same for all the dipoles. 
The correlation function is defined by 

where the Z^(x,y) is the partition function in the presence of the external charges, given by 

Zt(x,y) = ]T e -!^+€0")) J[ A(/2(j)). 

pL={PU),jeL} jeL 

where = £(6 jiX - 6 jiV ), and Z L = Z| =0 . 

Throughout this paper we suppose that the external charges are located far apart from each other, 
and that x and y do not belong to A. 

3 Electrostatic Inequality 

The first step to obtain an upper bound on the correlation function consists in driving down the 
activity of the dipoles, so as to obtain estimates valid for small values of z. This is done by the 
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electrostatic inequality. To apply this inequality, it is more convenient to define a dipole charge 
density by 

j,oi 

and write the electrostatic energy as E(p,i) = E(p) = (p, — A _1 p). If we carry out the summation 
on the dipole moment variables in (|2.1|) , we get 

P h& 

where z a (j) = ^(j)^ ^'. Because of the discrete orientation of the dipoles, p a (j) takes the values 
{0, ±1} in such a way that Pi(j)p2(j) = 0. 

Now we replace the charge density p, whose support lies in A, by a charge density p, whose support 
lies in A, in such a way that the interaction energy between the two renormalized densities remain 
unchanged, but that the self-energy of p is smaller than the self-energy of p. We use the difference 
in self-energies to renormalize the activities IQ. 
The renormalized charge density is given by 

P = p + W\' (3 ' 1} 

For the purpose of renormalizing the activity of the dipoles, we write 

Z = ^e-^e-^- E ^l[z Q (j). (3.2) 

P j,a 

The difference in self-energy is estimated as follows: 
So that we obtain 

Z = J2e-l E ^l[z a (j), (3.3) 

P j,a 

with the estimate z a (j) < z a (j) exp{-(/3/||A||)p^(j)}. 

Since the activities are all equal, all the renormalized activities can be bounded by the same constant 
(take 1 1 A|| = 8) : 

z a {j) <z = zexp{-^}, (3.4) 
which can be made arbitrarily small for f3 sufficiently large. 



4 The External Charges 



We apply the above reasoning to the partition function of the system in the presence of external 
charges, and obtain 

P j,a 

However, in this case we will need to renormalize the external charges as well. This is done by the 
shift 

a = p + Z+j(Aa)(j) 
where 7 depends on (3 in a way to be specified later, and 

a(j) = C(j,x)-C(j,y). (4.1) 

The function a(j) satisfies, for large \x — y\ and for \j — j'\ = 1, the bounds [3, chapter 7]: 

\a(j) - a(f)\ < const ( \ + 1 ) , (4.2) 

\\j-x\ + l \j-y\ + lj 

\a(j)-a(j')\ <const^J^, (4.3) 
From the definition of the Laplacian operator and (|4.1|) , we also have 

E Hj) -<f)f < (a,-Aa) =a(x) - a(y) w -\og\x-y\. (4.4) 
\j-j'\=i 71 

After the shift we get, following Eqs. ( |3~2| ) and ( |3.3| ), 

Z*(x,y) = e ( i-~ (){a '- Aa) Z 1 , (4.5) 

where 

^=E e " |i?(<j) n^o> 7Q(i) ' ( 4 - 6 ) 

with 

Ja(j) = lHaWoflij), (4-7) 

and the definition of a(j) follows from ( |3.1| ). 
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5 The Cluster Expansion 

To derive an upper bound on the partition function ( |4.6| ) we write 

n = n uj) cosh 7a (i) + r 

where (/x = 1, 2) 

R = 11(^0) cosh7 a (j) + z a (j) sinh7 Q (j)) - JJ z a (j) coshy a (j) (5.1) 
Thus ( |4.6| ) becomes 

= E e ~ f ^ I IT UJ) cosh 7a (j) + i? i (5.2) 

a { J> ) 

The function y a {j) depends on S a a(j). However, note that 

5 a a(j) = a(j) - a(j + e) = a(j) - a(j + e) when j x, y, 
since (Aa)(j) = for j ^ x,y. 

From the definition (|4.7|) , and the bounds (|4.2| ), ( |4.3| ) we see that the argument of cosh(-) is small 
for \j — x\ ^> 7 and \j — y\ ^> 7. Thus we can write 

cosh 7Q (j) - 1 < ~(1 + e) 7 2 [a(j) - a(/)f 

for some e > 0. 

For \j — x\ < 0(7), \j — y\ < 0(7), we estimate cosh(-) — 1 by a constant. 

In both cases we can take z sufficiently small so that the ( |5.1| ) term can be bounded by a constant 
for large j3. 

Under these circumstances, the i?-terms in ( |5.2|) can be neglected: 



Zi < KpJ2e~> E{(7) I[z a (j)cosh la (j) (5.3) 



for some positive constant Kp. 

Now we develop a cluster expansion for I\j a cosh7 Q (j) as follows 



I]!! + (cosh 7a (i) -l)] = l + S 1 + S 2 + S 1 S 2 (5.4) 



where 



with 



m a >l jl,...,jm a 



D a {ji) = cosh j a (ji) - 1. 



(5.5) 



(5.6) 



The sum in ( |5.5|) above runs over all subsequences of {ji, ...j ma }, m > 1. 
We can distribute the activities inside S a , to get 



Z x <K p \z + Y^ e~^\ Sl + s 2 + Sl s 2 ) | , 



(5.7) 



where s = Si + s 2 + Sis 2 and, as in (K 



S <*= J2 Za(jl)A*0'l)- • ■ Z a(jm a )D a (j ma ] 

m a >l jl,...,jm a 



Thus we see that Z\ admits the upper bound 



Z\ < KaZil + Si + S 2 + SiS 2 )- 



(5.J 



Here the s terms in (|5.7|) have been taken with fi a {j) = 1 for all j, a, and we have used that Z > 1. 
If we undo the cluster expansion (5^) we obtain 



Z 1 <K f ,Z]ftl + z a (j)D a (j)), 



(5.9) 



with fJ, a (j) — 1 f° r all j> a - 

Substituting (|5.9|) in ([4.6|) and using (|3.4j) and (f4.4|), we get 



The optimal choice for 7 is 7 = (3(1 + (1 + e)z/?)) _1 , for which we finally get 



(5.10) 



A. 1 



which concludes the argument. 
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